On the nonextensivity of the long range X-Y model 

Raul Tora|] 

Instituto Mediterrdneo de Estudios Avanzados (IMEDEA), 
CSIC-Universitat de les Illes Balears 
Ed. Mateu Orfila, Campus UIB, 07122 Palma de Mallorca, S'pcwfl 

(Dated: February 6, 2008) 

Abstract 

It will be given analytical and numerical evidence supporting that the X-Y model yields an 
extensive, i.e. proportional to the number of degrees of freedom N, internal energy U for any value 
of the interaction range. 
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The X-Y is a well known model of Statistical Mechanics. It is defined by a set of con- 
jugate pairs of variables (^i, i?«)i=i,...,jvj N being the number of degrees of freedom, with a 
Hamiltonian function 

H = \±* + Z± 1 -™f- 9 >^ K + V, (1) 

Z i=l Z i,i=l ' ij 

where the kinetic, K, and potential, V, contributions to the Hamiltonian have been iden- 
tified. In this expression is the distance between the sites of a regular lattice of 
dimensionality d and lattice spacing equal to one. For example, in the linear case that will 
be analysed in detail later, d — 1, it is = \i — j\. Periodic boundary conditions are 
assumed, and it is understood the convention of minimal distance between sites i and j. 
The constant k sets the interaction strength. The definition is such that, for ferromagnetic 
coupling k > (the only one considered in this paper) the minimum value of the energy is 
equal to zero. 

The parameter a sets the interaction range and its variation allows to recover some limits 
of interest. For example, the case a = oo, for which r~j a = unless = 1, is the nearest 
neighbor version, whereas the opposite limit, a = 0, is a mean field Hamiltonian in which 
all sites interact equally. This mean field limit is usually studied under the condition that 
the interaction strength is of order iV _1 , i.e. by setting k = e/N with e = 0(1). It is known 
that in this case the model exhibits a phase transition from an ordered to a disordered phase 
at a critical value (3 C = 2/e. This result follows 2] from an analysis of the partition function 

/oo roc r2iT r2n 

d Pl ... dp N <f0i.../ d6 N e- im (2) 
-oo J-oo JO JO 

at inverse temperature (3 = l/fc^T, as usual. 

The model has received recently some attention 3} for a value of the interaction strength 
k = 0(1) as an example of a non-extensive system. It has been argued^] that in this case 
the internal energy U = (T~l), as well as other thermodynamical potentials, should scale as 
U ~ N 2 ~ a / d for a < d (the so-called nonextensive regime), whereas for a > d, the extensive 
regime U ~ iV is recovered. 

In this paper, I will present arguments showing that the internal energy is an extensive 
function U (N, T) = Nu(T) for any value of a. The result is based upon an exact calculation 
of the partition function in the mean-field limit a = as well as some numerical simulations. 
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Let me first sketch the proof for a = 0. Introducing the total "spin vector" M = J2i 7 ^i 
with rhi = (cos 0i, sin 6i) it is possible to write the potential energy as V = |(iV 2 — M 2 ). Us- 
ing this relation and the Hubbard- Stratonovich transformation, one arrives at the following 
expression for the partition function Q]: 

JV 

z-ffi^QgL^+MSr-*** (3) 

where Zk = yf) 2 an d Zy are, respectively, the kinetic and potential contributions and 
I n (y) is the modified Bessel function of order n. The standard calculation^] assumes that 
k = O(N) and finds the dominant saddle point contribution to the above integral in the 
limit N —>■ oo. A similar calculation can be carried out assuming that k = 0(1). In this 
case, the location of the saddle point yo is given as the solution of 

1 22/o , h(yc? 



+ T7^( = 0. (4) 



Ny (3kN I Q (y Q ) 

This can be solved in the limit iV — > oo, with the result that yo = — | + 0(N^ 1 ). The 
resulting partition function is: 

/2tt\ 2 1 / 2tt \ 2 

It turns out that the Helmholtz free energy F = — ksT In Z is not an extensive function 
since it contains terms of order iVlniV. However, these terms disappear when computing 

iln. 
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the internal energy U = — dl a B Z - This yields U = (N + l)//3. In the thermodynamic limit 



the dominant term is: 

showing that the internal energy is indeed an extensive function of the number of degrees of 

freedom. Note that the kinetic and potential terms in the Hamiltonian contribute equally, 
N 

(K) = (V) = — , to the internal energy. 
Zp 

The nearest neighbor interaction, equivalent to the limit a = oo, has been studied in 
Ref. 0, Q| for the linear lattice d = 1. It is found that the internal energy is again an 
extensive function, namely: 



U = N 



4+^1 



2(3 \ Iq(Pk) 



(7) 
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Therefore, since the limits a = and a = oo both yield an extensive internal energy, it is 
natural to speculate that a similar behavior will hold for any value of the parameter a. In 
order to check this result I have performed Monte Carlo simulations of the one dimensional 
system given in for the values of a = 0, 0.5, oo and a value of the coupling strength 
k = 1. Whereas the limiting cases a = 0, oo are only for the purposes of checking that the 
numerical simulations do reproduce the theoretical results, the intermediate value a = 0.5 
can not be compared with any a priori theoretical calculation. 

I have used a standard Metropolis algorithm j| to sample the phase space with a prob- 
ability proportional to e _/3W : (i) Choose randomly a lattice site i. (ii) Propose a change 
(pi,9i) — > (PijO'j) of the variables associated to this site; the proposed values for p[ and 6[ 
are selected randomly from a uniform distribution in (pi — 5 p ,pi + S p ) and ($i — Sg, 9i + 5g), 
respectively, (iii) Accept the change with a probability min[l,e~^ ], being AH. = 7i' — TH. 
the change in energy involved in the proposed change. The values of 6 P and 6g are chosen 
such that acceptance probability remains close to 50%. Typically, after 10 5 Monte Carlo 
steps per particle for equilibration, results have been averaged for 10 6 Monte Carlo steps per 
particle for each value of the inverse temperature (3. In all cases, the following system sizes 
have been taken: N = 50, 100, 200. 

Figures EEl plot the internal energy rescaled by the number of degrees of freedom, U/N, 
as a function of the inverse temperature j3, for the differer values of a considered.. Figure 0] 
plots the rescaled potential energy, V/N, in the case a = 0.5. In all cases the data for all 
values of iV collapse in a single curve, showing the extensive nature of the internal energy 
for all the values of a. Moreover, the cases a = and a = oo, Figs. [T] and follow 
the theoretical predictions given by Eqs. (0) and (jjj), respectively. Note that, as shown in 
Figs. El there is not much quantitative difference between the values of the internal energies 
corresponding to these two limiting cases. 

As other authors have found it is also possible to map the results of the X-Y model 
with k = 0(1) to the same model with k = O(N) if one allows to rescale the temperature 
T by the number of degrees of freedom N. In this way it is possible to achieve a sort of 
non-extensive regime in which U(N, T) = N 2 u(T/N) (for a = 0). The results of this paper, 
however, concern the true thermodynamic limit, iV — > oo, while keeping T = 0(1). It 
has been shown that this limit yields an extensive internal energy and that this extensive 
behavior can be recovered by simple Monte Carlo simulations. It is interesting to remark 
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that the same extensive behavior will be found in microcanonical simulations. In fact, it is 
possible to compute the density of states g(E) for the potential part of the Hamiltonian as 
the inverse Laplace transform of the partition function ^y[9|. This gives: 



where c is a constant of order one. As shown in figure |5J this result is fully consistent 
with a numerical calculation of the density of states 10| which predicted empirically a form 
g(E) = e N ^( E / N2 )_ Within this microcanonical ensemble, the internal (potential) energy is 
found by using (3 = ^Mpl or E = N/2/3, the same result we obtain before for the average 
potential energy (V). 

The existence of this extensive regime can be considered surprising at first sight, since the 
interaction terms are very strong. From the mathematical point of view, the extensive regime 
is due to the continuous nature of the excitations of the ground state of the Hamiltonian 
(JIJ. A similar analysis for the long range Ising model in which the spins can only vary by 
an integer amount indeed produces a non-extensive internal energy. In the X-Y model with 
long range interactions and k = 0(1) there is no phase transition at any finite value of the 
temperature in the thermodynamic limit, and the system is always in an ordered situation. 
In fact, the order parameter is M = O(N) for all values of the temperature. 

In summary, it has been given analytical and numerical evidence showing that the internal 
energy derived from the X-Y Hamiltonian is an extensive function of the number of degrees 
of freedom, U(N,T) = Nu(T), for any value of the parameter a setting the interaction 
range. Other thermodynamic potentials, such as the Helmholtz free energy, contain terms 
of the order N log N. Non extensive regimes require to scale the temperature with N. 
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FIG. 1: Internal energy U rescaled by the number of degrees of freedom for a = and coupling 
strength k = 1. The symbols are the result of the Monte Carlo simulations, whereas the solid line 
is the theoretical prediction Eq. ([5]). Circles, diamonds and triangles correspond, respectively to 
N = 50, 100, 200. It is observed a very good collapse of the data for all values of N showing the 
extensive character of the internal energy. 
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FIG. 2: Similar as Fig. ^for a = oo. The symbols are the results of the numerical simulations 
(same symbol meaning than in FigQ) and the solid line is the theoretical prediction Eq.([7jl. 




FIG. 3: Similar as Fig. ^for a = 0.5. The symbols are the results of the numerical simulations 
(same symbol meaning than in Fig^). For comparison, the theoretical predictions for a = 0, Eq. 
iJSJ), solid line, and a = oo, Eq.Q, dotted line, are also included. 
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FIG. 4: The potential part V of the internal energy in the same case than in figure |21 The symbols 
are the results of the numerical simulations (same symbol meaning than in Fig^) whereas the solid 
and dotted lines are the theoretical predictions for the limit cases a = and a = oo, respectively. 




FIG. 5: Density of states g(E) for the potential part of the X-Y Hamiltonian in the case a = 0. 
The solid line is the theoretical prediction Eq. @, with a constant c = 2.54, while the dots are 
the result of a numerical calculation ^(| of the same quantity. Please note that in reference ^(| 
the data were plotted using base 10 logarithms instead of the natural logarithms used here. 
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